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Abstract 

We define the collinear factorization scheme, which absorbs only the collinear physics into the 
parton distribution functions. In order to isolate the collinear physics, we introduce a procedure 
to combine real and virtual corrections, canceling infrared singularities prior to integration. In the 
collinear scheme, the factorization scale \x has a simple physical interpretation as a collinear cutoff. 
We present a method for choosing the factorization scale and apply it to the Drell-Yan process; 
we find fi ~ Q/2, where Q is the vector-boson invariant mass. We show that, for a wide variety of 
collision energies and Q, the radiative corrections are small in the collinear scheme for this choice 
of factorization scale. 
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I. INTRODUCTION 



In the parton model, a hadron is regarded as a collection of quarks, antiquarks, and 
gluons, each of which carries some fraction x of the hadron's momentum, with a number 
density f{x, //), where fi is the factorization scale. Qualitatively, the factorization scale 
corresponds to the resolution with which the hadron is being probed. To calculate the cross 
section for processes in hadron-hadron or lepton-hadron collisions, the partonic cross section 
is convolved with the corresponding parton distribution functions f(x,fi). The partonic 
(hard-scattering) cross section is independent of the factorization scale fi at leading order in 
perturbative QCD, but depends logarithmically on /i at next-to-leading order and higher. 
When calculated to all orders in perturbative QCD, the hadronic cross section is independent 
of \i. However, at any finite order in perturbation theory, the calculated hadronic cross 
section depends on /z. This dependence is usually significant at low orders in perturbation 
theory. One could conclude that the only way to obtain a reliable prediction is to calculate 
higher-order corrections until the factorization-scale dependence is reduced. However, there 
are many processes for which higher-order corrections are not available. It is in these cases 
that the question of factorization-scale choice is most important. 

If the hard-scattering cross section is characterized by a single scale Q (such as the 
invariant mass of the lepton pair in Drell-Yan production), then the factorization scale 
is usually chosen to be of order Q, simply because there is no other scale in the problem. 
However, this reveals only the order of magnitude of /i. It is common to vary the factorization 
scale over some interval, perhaps Q/2 < fi < 2Q, but there is no objective argument for 
either the central value or the range of variation. If the hard-scattering cross section depends 
on more than one scale, then the choice of factorization scale becomes an even murkier issue. 

One argument against trying to do better than simply choosing // ~ Q (we use ~ to 
denote order-of-magnitude equality) goes as follows. Since the hadronic cross section is a 
physical quantity that does not depend on any factorization scale, the factorization scale 
is unphysical, and therefore one cannot make a physical argument for its choice. It is in 
part due to arguments such as this that there has been little effort to try to do better than 
choosing /x ~ Q. 1 

1 One attempt is based on the principle of minimal sensitivity [l| , and another on complete renormalization 
group improvement Q . 
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On the other hand, the statement is often made that the factorization scale separates the 
short- distance physics of the hard-scattering cross section from the long-distance hadronic 
physics [3]. The qualitative statement made earlier, that \i corresponds to the resolution 
with which the hadron is probed, falls into this class. This is contrary to the attitude that 
the factorization scale is unphysical. 

In this paper we put forward a physical argument for the choice of factorization scale. 
We restrict our attention to Drell-Yan production, but the argument can and hopefully will 
be extended to other processes in the future. 

The physical argument we advance is that the factorization scale should be chosen such 
that collinear (long-distance) physics is included in the parton distribution functions, and 
noncollinear (short-distance) physics in the hard-scattering cross section. This is not a new 
idea, and goes back to the origins of the parton model. 2 What is new is the implementation 
of this idea in practice. We introduce a new factorization scheme, which we dub the collinear 
scheme, in which the factorization scale corresponds to a collinear cutoff. We then argue for 
a method to choose the factorization scale in that scheme. 



The method we pursue here was 



irst proposed in the context of Higgs-boson production 



in association with bottom quarks 5), Q] , and was refined and elaborated upon in Ref . [?J ■ 
An alternative approach in the same spirit has been developed in Ref. 8y. 

All of the above studies are in the context of Higgs-boson (both charged and neutral) 
production via a bottom-quark distribution function, which is a rather exotic process, both 
in the initial and the final state. However, the ideas developed there are of general validity, 
and should be applicable to all parton-model calculations. In this paper we further develop 
7j and apply it to one of the most basic processes of perturbative QCD, 



the method of Ref. 



namely Drell-Yan production 



101 ] . Our hope is that this method can be generalized 
to all parton-model calculations, finally satisfying the desire for a systematic method of 
factorization-scale choice. 



See, for example, Ref. 
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II. THE COLLINEAR FACTORIZATION SCHEME 



The parton distribution functions are evolved from one factorization scale to another 
via the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. These equations 
sum collinear logarithms into the parton distribution functions. In this section we define a 
factorization scheme in which the factorization scale \x has the interpretation of a cutoff in 
the integration over the virtuality of a propagator associated with collinear radiation. 




V 



FIG. 1: Drell-Yan production of an electroweak gauge boson V. 

Consider the Drell-Yan process, shown in Fig. [TJ The colliding quark and antiquark q~j 
annihilate into an electroweak gauge boson V of invariant mass Q; we do not consider the 
subsequent decay of the boson into lepton pairs, which is irrelevant to our discussion. The 
leading-order (LO) cross section is 

Co) 47r 2 a _ 

= c iM^ ® Qoj + qoj ® qoi){zo), (i) 

where VS is the hadronic center-of-mass energy, z = Q 2 /S, (shown in Tabled]) specifies 
the coupling of the quarks to the boson V, and the convolution fx ® f 2 of two functions fi 
and f 2 is defined by 

(/i ® h){x) = I dxx dx 2 S(xix 2 - x)/i(xi)/ 2 (x 2 ). (2) 
Jo Jo 

We will also need this cross section calculated in dimensional regularization with 4 — 2e 
spacetime dimensions (/id is the 't Hooft mass): 

Oft = ~ e )^n Yl C iMoi ® Qoj + Qoj ® 9w)(«b)- (3) 
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TABLE I: Couplings of the quarks to electroweak gauge bosons. Here syy = svaOw, c\v = cosOw, 
and Vij denotes an element of the Cabibbo-Kobayashi-Maskawa matrix. 
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A. Initial gluons 




FIG. 2: Correction to the production of an electroweak gauge boson V due to initial gluons. 

The quark distribution function receives corrections due to the splitting of a parton 
into a pair of collinear partons, one of which is a quark. In order to find the correction 
to qi(x) arising from the splitting g — > qiqi, we consider the next-to-leading-order (NLO) 
correction to the Drell-Yan cross section due to the process gqj — > q~iV, shown in Fig. [2j 

Let s = (pi +P2) 2 , t = (pi —P3) 2 , and u = (j>2 —pz) 2 be the usual Mandelstam variables, 
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with s + t + u = Q 2 . The cross section for this process is 

-S Ar , rO rO 
'Q 2 



(i) iraa. 

(7 - = 

m 3S 



rS j rO rO 

VCy — / dt du5{s + t + u-Q 2 ) 

"j JQ 2 s J-oo J-oo 



x (9 ® Qj + Qj ® #) 



s -t 2uQ 2 



(4) 



-t s st 

This expression manifestly shows the singular behavior near t = 0, which corresponds to the 
splitting of the initial gluon into a collinear quark-antiquark pair. The integrand contains a 
simple pole at t — 0; thus this expression for the cross section is divergent. 

If we use the delta function to perform the integration over u, we can break the cross 
section into two pieces: a "collinear" piece, comprising the terms proportional to 1/t, and a 
"noncollinear" piece, comprising all other terms. Thus = (Tg^°°^ + Ggq > where 

(i,coi) _ 27rao; s ^ n f S ds f° dt r„™* - 0n ^s^ D ( Q 2 

O n 



m 3S 

hi 



E Q ' I t ^ (9 8 * + * 9) (I) p - (?) ' (5) 



(i,non) Traa, X^ri f ds f mi - , - ^ ^ s A 2< 2 ~ * /«\ 
^ -^-Z^ C v — dt{g®qj + qj®g)[-\ , (6) 

and Pqg(z) = \\z 2 + (1 — z) 2 ] is the DGLAP splitting function. The differential cross 
sections —tda^' col ^/dt and —tda^ /dt corresponding to these two pieces (combined with 
their counterparts corresponding to the analogous process qig — > Vq{) are shown, for the case 
of real Z-boson production at the Tevatron, in Fig. [3j 3 The two curves resemble a "plateau" 
and a "hump," respectively, and their contributions to the cross section are proportional 
to the areas under the curves. The noncollinear contribution is small, while the collinear 
divergence is evident in the plateau, which extends infinitely far to the left. We must absorb 
this divergence into the parton distribution functions if we are to obtain a finite prediction 
for the cross section. 

In order to do this, we define the collinear factorization scheme as follows. We define 
a counterterm a g l q by imposing a cutoff — fi 2 on the integration over t and replacing the 



3 We take the renormalization scale for a s to be Q, here and throughout this paper. This is a separate 
issue from the choice of the factorization scale. 
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FIG. 3: The quantities —tda g ^ q g/dt and —tda g ^°g/dt, defined via Eqs. © and ©, for the case 
of real Z-boson production at the Tevatron. The "collinear" curve, —tda^^g/dt, passes through 
50% of its limiting value when = 0.53Q. 

integrand with 1/t times the residue of the collinear pole: 




The counterterm is determined solely by the collinear part of the cross section; Eq. ([8]) is 
identical to Eq. (j^D except for the replacement s—Q 2 — > fx 2 in the lower limit of the t-integral. 
We regard /i as the factorization scale in this scheme. By subtracting the counterterm from 
the bare cross section, we obtain the explicit correction in the collinear scheme: 
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By a rearrangement of terms, 



99 91 



2naa s 
3S 



hi 



ds ( s 

— {g® qj + q j ®g)(- d 

2 S \ O 



X 



dt 



s+Q 2 —t 



— P qg [ — ) + 



dt 



2Q 2 



s+Q 2 



2s 2 



(10) 



we obtain an expression in which the integrals are finite. The result is 



a 



a 



<j<i 



27raa s 
3S 



J2 C v j ~ & ® % + Qj ® 9) ( 



P w (z)ln^^ + i(l + 2z-3z 2 ) 



If an explicit expression for the counterterm is desired, it is necessary to regulate 
the collinear divergence. This can be done by working in dimensional regularization and 
suitably generalizing the definition of the counterterm: 

«o 



4> = 



lim 



-t) 



da {i y 

i+t uu gq 



dt 



dt 



-a 2 (~t) 



1+e 



(12) 



We perform the integration over t and expand in powers of e to obtain 



a 



(i) 

91 



2-KOLOls 

3S 



— {g® Qj + qj ®g)( 

■° z 



X 



47T/i 2 



P qg (z) - z(l - z) 



(13) 



This result can be expressed concisely as 

47r 2 a 



4> 



3S 



(14) 



where 



<%(x) = - 



2tt 



1 cfz /x 

— 5 (- 

z \z 



1 , An/A 



P 3ff (z) - z(l - z) 



(15) 



The correction 8q~j(x) to the antiquark distribution function, due to the process — > Vqj, 
is identical. These corrections are absorbed into the parton distribution functions if we 
define qi{x) = qoi(x) +5qi(x) and q~j(x) = q~oj(x) +dq~j(x), and calculate the LO cross section, 
Eq. d2J), with qi(x) and q~j(x) in place of qoi(x) and q~oj(x). The explicit correction, a^j — 0gq, 
calculated in dimensional regularization, matches Eq. pip . 
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We come now to the question of the choice of /x. It is possible to choose /x such that 



(i) 

91 



(l.col) 
9<? 



With this choice, the collinear physics is entirely absorbed into the parton 



distribution functions, and the explicit correction to the cross section is equal to er 



(l,non) 

gq 



However, a precise determination of the value of /x that accomplishes this is unnecessary; in 
practice one can come close by choosing /x to equal the value of y/—i for which —t da^ ' /dt 
passes through 50% of its limiting value. This is the convention we will follow in this paper. 
For the case of real Z-boson production at the Tevatron (Fig. [3]), this prescription indicates 
a scale \i = 0.53Q. 
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FIG. 4: The factorization-scale dependence in the collinear scheme of the LO and initial- gluon 
NLO contributions to real Z-boson production at the Tevatron. The factorization scale indicated 
by the plateau is 0.53Q. 

Fig. H] shows the factorization-scale dependence of the explicit correction to the cross 
section in the collinear scheme, given by Eq. fTTTl) (together with the analogous contribution 
from q^g — > Vqj), for the case of real Z-boson production at the Tevatron. This correction 
is small near the scale fi = 0.53Q, which supports our argument for this scale. However, 
the correction shown in Fig. H]is also small for /x = Q, so we have not yet demonstrated the 
superiority of using a scale other than Q. This will become evident when we consider higher 
values of Q in Section IIHI 
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B. Real and virtual gluons 



Pi 




P1-P3 




+ 




FIG. 5: Correction to the production of an electroweak gauge boson V due to real gluon emission. 

Now consider the NLO correction due to gluon radiation, shown in Fig. [5j The cross 
section is 



(l.real) _ 

w " = 9S 



dt / w ■ >-< 

oo J — oo 



+ t + u - Q z 

x (ft ® ?j + ?j <g> ft) ^) 



't it 2sQ 2 
_ + - + — ^ 
u t tu 



Using the identity 
we rewrite this as 

(l,real) _ S7iaa s 



1 



1 



11 



95 ^ lJ 



1 1 

tw s — Q 2 ^— t — u , 



30 (is /"° Z" 
— / dt I du5(s + t + u-Q 2 ) 

2 S J —oo J -oo 



(16) 



(17) 



x (ft ft + ft <g> ft) 



s 2 + g 4 / 1 
s-g 2 v-t 



— 1-2 

— u 



;is) 



There are two collinear singularities, at £ = and at it = 0, corresponding to the emission 
of a collinear gluon by the initial quark and antiquark, respectively. Each is accompanied 
by an infrared divergence l/(s — Q 2 ). This soft singularity resides at the intersection of 
the two collinear singularities (t = u = 0) and threatens to obscure the collinear physics in 
which we are interested. It is cancelled, however, by the infrared divergence from diagrams 
involving virtual gluons. We therefore turn our attention to these, beginning with the vertex 
correction shown in Fig. [6j 
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FIG. 6: Vertex correction to the production of an electroweak gauge boson V. 

The inclusion of this diagram modifies the tree-level amplitude by the replacement 7 M 
^[l + 5F 1 (Q 2 )], where 



16 



d A k 



k 2 (k - Pl ) 2 (k + P2 y (k 2 - A 2 )(k - Pl ) 2 (k + P2 y 



9) 



In the second term we have introduced the Pauli-Villars regulator A >Q. The loop integra- 
tion is simplified by introducing Feynman parameters x, y, z and shifting the loop momentum 
to £ = k — xpi + yp2- The result, up to terms that vanish in the A — > oo limit, is [111 ] 

32 



x 



-ma. 



d 4 



/ dx dy dz 5(x + y + z — 1 
Jo Jo Jo 

2 -2{z + xy)Q 2 I 2 



(20) 



(2^- 



(£ 2 + xyQ 2 ) 



(£ 2 - zA 2 )\ 



2a, 
3tt 



/ dx l dy I dz 8{x + y + z — 1) 
Jo Jo Jo 



(21) 



x 



1 



1 - In 



zA 2 



1 — z \x ' y J \ xyQ 2 

The integrand contains terms proportional to 1/x and 1/y, which correspond to collinear 
singularities at x = and at y = 0. Each is accompanied by a factor 1/(1 — z); this is 
the infrared singularity which cancels against that of the real-gluon contribution. The term 
involving In A is an ultraviolet divergence. 

We must also include the correction due to wavefunction renormalization of the external 
quark lines, shown in Fig. [71 This involves the quark self-energy E(p), which for mass- 
less quarks takes the form S(p) = T,'(p 2 )p. The tree-level amplitude is modified by the 



multiplicative factor 1/yl + E'(0) for each external leg. 
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V + 




FIG. 7: Wavefunction renormalization correction to the production of an electroweak gauge boson 
V. 

Let us first consider the correction to the incoming quark. The self-energy is 



£(Pi) = yTT^s J 



d A k 



k 2 (k - pi) 2 (k 2 - A 2 )(k - Pl ) 



(22) 



This integral contains an ultraviolet divergence (regulated by the inclusion of the second 
term), which cancels against that of the vertex correction. It also contains a collinear 
divergence, which we want to express in a form that can be easily combined with those we 
have encountered in the vertex correction. The natural way to proceed with the calculation 
of is to introduce Feynman parameters x and y, with x + y — 1. This will not suffice 

for our present purposes; we need a two-dimensional Feynman-parameter space in order to 
match the vertex correction. We can achieve this by multiplying by unity: 

^-f l ){k + p 2 ) 2 ^-p i )(k + p 2 ) 2 



S(pi) 



32 



-irta q 



(Pk 



k 2 {k-pi) 2 {k+p2) 2 (k 2 - K 2 )(k - pi) 2 {k + p 2 f 



(23) 



The denominator now matches the denominator of the vertex correction, and we can there- 
fore introduce matching Feynman parameters. 

We can think of p 2 , which has appeared out of left field, as the momentum of the incoming 
antiquark in q^qj — > V. We take the antiquark and the vector boson to be on shell, i.e. p\ = 
and (p\ + p 2 ) 2 = Q 2 , but p\ must be kept off shell for the moment. 

We introduce Feynman parameters x, y, z and shift the loop momentum to t = k — xpi + 
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yp2- 



MPl) = ^Kia s 



dx I dy 



x 



dzS(x + y + z — 1) 

1 1 

[£ 2 + xyQ 2 + xzp\f ~ (£ 2 -zk 2 f 



x | [(2 - 3x)f l - (1 - 3y)^]£ 2 + 2x[(l - x)^ + y#J [(1 - y)Q 2 - zp\]\. (24) 

The presence of terms proportional to fl 2 may be surprising, but these terms integrate to 
zero and do not contain a collinear singularity. Thus they need trouble us no further. At 
this point we can read off the coefficient of p and set p\ = 0: 

£'(())= — —Tti<y s f dx f dy f dz5(x + y + z — l) 
<J Jo Jo Jo 

'(2 - 3x)£ 2 + 2x(l - x)(l - y)Q 2 (2 - 3x)£ 2 



(25) 
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2a 



dH 



(2vr) 4 
• l 



zA 2 
xyQ' 



(26) 



(£ 2 + xyQ 2 ) 3 (£ 2 -zA 2 )\ 

/ dx dy dz 5(x + y + z — 1) — h x + (2 — 3x) In 
3vr Jo Jo Jo [y 

The term involving In A is the expected ultraviolet divergence. There is also a single collinear 
singularity at y — 0. 

The self-energy correction to the incoming antiquark line is completely analogous, but 
it is appropriate to interchange the parameters x and y. The order-o; s contribution to the 
QiQj ~^ V amplitude due to virtual gluons is therefore equal to the tree-level amplitude times 
the factor 



a = ^(q 2 ) - ^£'(0) - itf(o) 



— — dx dy dz 5(x + y + z — 1) 
3tt Jo Jo Jo 

z(3 - z 



(27) 
(28) 



- + -)+3-z-(l-3z)ln 

re y 



zA 2 



This expression still contains the ultraviolet-divergent term ln(A 2 /Q 2 ). However, its coef- 
ficient integrates to zero, so that we can discard this term. The NLO contribution to the 
cross section is then 



r (l,virt) 



■qq 



(2 Re A) a, 

8naa 



(0) 

gq 



(29) 
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9S 

z{3 - 



Cij(qi ® qj + qj <g> gj)(2 ) / dz / dy / d^ 5(x + y + z - 1) (30) 

Jo Jo Jo 



1 1, 

-+- ) +3 

x y / 



[l-3z) In — 
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This must be combined with the correction due to real-gluon emission to yield an infrared- 
finite result. 

The hadronic real-gluon-emission cross section involves an integration over the Mandel- 
stam variables s, t, and u, subject to the constraint s + t + u = Q 2 . There are collinear 
singularities at t = and at u — 0, and where the two meet there is an infrared singularity at 
s = Q 2 . Similarly, the virtual-gluon cross section involves an integration over the Feynman 
parameters x, y, and z, constrained by x + y + z — 1. There are collinear singularities at 
x = and at y = 0, and where the two meet there is an infrared singularity at z = 1. The 
two spaces can be mapped onto one another in such a way that the singularity structures 
match exactly: 




FIG. 8: The region of Feynman-parameter space bounded by x = 0, y = 0, and z = is mapped 
by Eq. (|31|) onto the region of Mandelstam-variable space bounded by s = oo, t = 0, and u = 0. 
The real-gluon contribution is bounded by s = S (or z = zq), while the virtual-gluon contribution 
extends all the way to z = (or s = oo). The heavy lines represent the collinear singularities, and 
the black dots represent the infrared singularity. 

We thus change variables from x, y, z to s, t, u, using the relation 

f dx [ dy [ dz5(x + y + z-l) =Q 2 [ ^ / dt f du5(s + t + u - Q 2 ). (32) 

JO JO Jo Jq 2 s J -oo J-oo 
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The virtual-gluon cross section becomes 

M^)-- 87[a ^J2C lJ (q^q 3 + q 3 ^q t )(z ) H% f dt f du 8{s + t + u - Q 2 

Jo 2 s J-oo J-oo 
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9S 



Q 2 



Q 2 (3s-Q 2 ) / 1 



1 \ 3s -Q 2 s-3Q\ sQ 2 

— + In — - 

—11 I s s tu 



(33) 



s-Q 2 \-t 

In this equation the symbols s, t, and u, while they have been designed to resemble Man- 
delstam variables, are not defined in terms of external momenta, but rather as particular 
combinations of Feynman parameters given by Eq. (1511) . In particular s ^ (p± + P2) 2 = Q 2 ■ 
We add Eqs. ffTBT) and (133]) to obtain the total cross section due to real and virtual gluons: 



(i) 
11 



9S 



'■J 



ds 



2 



dt 



du S(s + t + u- Q 2 



-(qi®qj + qj®q: 



x y(S - s)(qi ® qj + q 3 ® q % ) {^j 
Q 2 \ Q 2 \Q 2 (3s-Q 2 ] 



s-Q 2 \-t -u 



S 



s-Q 2 



1 1 \ 3s - Q 2 s-3Q 2 sQ 21 

— + + In — 

—t —u Is s tu 



(34) 



The infrared cancellation is now manifest, as the residue of the s = Q 2 pole vanishes. The 
cancellation of the infrared divergence prior to integration is reminiscent of the work of 
Refs. |l2,|3,Q|. 

As in Eqs. (JSJ) and (jSJ), we decompose into collinear and noncollinear pieces: = 



a 



(l,col,t) 



11 



(l,col,u) 

11 



17 aa " - o-S' n ° n) , where 



(l,col,t) 
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inaa s 
9S 



ds 



2 s 



dt 



9(S - s)(qi <g> qj + qj ® q, 
- {q% ® qj + q~j ® qi) 



s\ s 2 + Q 4 



SJ s(s-Q 2 
Q 2 \ Q\3s-Q 2 ) 



S J s 2 (s - Q 2 



(35) 
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(l,non) 
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8iraa t 
9S 



^Qj [ ^ f dt \ du5{s + t + u-Q 2 ) 

~~ JQ 2 S J-oo J-oo 



xh(S- s)( fi ® qj + qj ® qi ) (|) [-2] 



- (qi ® qj + qj ® qi) [ -j- J — 



Q 2 \ Q 



3s -Q 2 s-3Q 2 sQ 



In 



tu 



(36) 



and ail' col ' u ^ is defined similarly to er£' c . The differential cross sections —tda& col '^/dt 
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and —tdaZq^ /dt are shown, for the case of real Z-boson production at the Tevatron, in 
Fig. E 
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FIG. 9: The quantities —t da qq -' co ^ / dt and —tda^ mm '/dt, defined via Eqs. ([35]) and ([36]) . for the 
case of real Z-boson production at the Tevatron. The "collinear" curve, —t da^' c ° l ' t ' / dt, passes 
through 50% of its limiting value when \f—i = 0.57Q. 



The t-channel collinear- scheme counterterm takes the form 
^ f°° ds f° dt 

2>« L 2 T / 



- (ft ® qj + qj ® ft) 



2naa s 
3S 



S J s 2 (s-Q 2 ) 
£ f £ ( ft s g, + g, ® ft ) (|) P g? (£) f ^, 



(37) 



(38) 



where P gg (^) = |[(1 + z 2 )/(l 



is the DGLAP splitting function. The w-channel 



counterterm cx^'"' 1 is similar. By subtracting these counterterms before integrating over t 
and u, as in Eqs. (Q and (jTUl) . we arrive at the following result for the explicit correction in 
the collinear scheme: 



(1) _(l,col,t) -(l.col.u) 
w qq qq qq 



2naa s 
3S 



'■J 
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2P qq {z) hi^ r + -(l + z i 



//- 



V 1-2 



- + 3j 5(1 ^ } 
81 + ^ 2 8 ■ 

3T^7 ln "-3 (1 - z) 



(39) 
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An explicit expression for the counterterms can be obtained by working in dimensional 
regularization: 

27raa s N „,v^ „ C 1 dz 



- fi,t) 

<T - 



_(l,u) 

(T - 

u qq 



3S 
1 



1 - e Vi? X] ^ / — (ft ® + * ® f ~ 

ln(l - * 



- 7 + In 



'■j 



1 + z^ 



4tt 2 



(40) 



There is one added complication in dimensional regularization: it is necessary to correct for 
double counting in the region where both t and u are near zero (see Fig. [TP]) . Therefore, 
after subtracting the two counterterms, we must add the quantity 



*W = 

u qq — 



lim 

t,u->0 








du ■ 



1 dz 



J2 C vf + (7) ^y 6 ^-*)- ( 41 ) 



By following this procedure, we recover the explicit correction given in Eq. (1391) 151 ] . The 
collinear contribution absorbed into the quark distribution function is 

47r 2 a 



!=(i) 

— (7 - 

2 qq 



35 



1 



X) ^ ^ ® * + * ® 5 ^ Cy) ' ( 42 ) 



where 



-7T- / — % 



1 , 47T/i 

- - 7 + In — - 

e /i^ 



4 
3 

14 \ 



Pqq(z)--(1 + Z 

8tt 2 



ln(l - z) 
1 - z 

5{l-z)--(l-z) 



(43) 



9 3 7" v " '' 3' 
The correction 5q~j(x) to the antiquark distribution function is analogous. 

The 50% rule described in Section III applied to the plateau curve in Fig. [9], indicates a 
factorization scale \x = 0.57Q, about the same as that indicated by the initial-gluon plateau. 
Fig. E] shows the factorization-scale dependence of both NLO corrections to the Drell-Yan 
cross section in the collinear scheme. At /i ~ Q/2, the correction due to real and virtual 
gluons is small, about —7% of the LO cross section. This correction is even smaller at 
/i ~ Q, essentially vanishing. However, there is no reason to expect the correction to vanish 
at the appropriate factorization scale; it is enough that it is small, indicating a convergent 
perturbation series. 
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FIG. 10: The counterterms corresponding to the t- and u-channel collinear singularities are formed 
by integrating over the shaded strips. We must be careful to avoid double-counting in the region 
where the two strips overlap. 

The collinear-scheme counterterms for the Drell-Yan process can be expressed simply in 
terms of convolutions of parton distribution functions with DGLAP splitting functions, as 
in Eqs. (jHJ) and ( 1381) . This suggests that the counterterms for other processes will have a 
similar form, and that NLO parton distribution functions denned in the collinear scheme 
will be universal. 

III. RESULTS 

In this section we show the factorization-scale dependence of the LO and NLO contribu- 
tions to the Drell-Yan cross section, in the collinear scheme, for several values of Q and yS. 
In some cases the initial-gluon correction has been multiplied by a large factor in order to 
make it more visible. 

Fig. [12] shows the case of virtual Z-boson production at the CERN Large Hadron Collider 
(LHC) (y/S = 14 TeV pp), with Q = 650 GeV. This value of Q has been chosen so that 
the ratio Q/y/S is the same as for the case of real Z-boson production at the Tevatron, 
presented in Fig[TTJ The factorization scales indicated by the plateaux are about the same 
as in that case, still in the vicinity of n ~ Q/2. Again, both NLO corrections are small near 
this scale. 
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FIG. 11: The factorization-scale dependence in the collinear scheme of the LO and NLO contri- 
butions to real Z-boson production at the Tevatron. The factorization scales indicated by the 
plateaux are 0.53(3 f° r initial gluons and 0.57(3 for real and virtual gluons. 
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FIG. 12: The factorization-scale dependence in the collinear scheme of the LO and NLO contri- 
butions to virtual Z-boson production at the LHC, with Q = 650 GeV. The factorization scales 
indicated by the plateaux are 0.52(3 for initial gluons and 0.61(3 for real and virtual gluons. 
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FIG. 13: The factorization-scale dependence in the collinear scheme of the LO and NLO contribu- 
tions to virtual Z-boson production at the Tevatron, with Q = 490 GeV. The factorization scales 
indicated by the plateaux are 0.37Q for initial gluons and 0.63Q for real and virtual gluons. 
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FIG. 14: The factorization-scale dependence in the collinear scheme of the LO and NLO contri- 
butions to virtual Z-boson production at the LHC, with Q = 3.5 TeV. The factorization scales 
indicated by the plateaux are 0.37Q for initial gluons and 0.58(3 f° r rea l an d virtual gluons. 
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FIG. 15: The factorization-scale dependence in the collinear scheme of the LO and NLO contribu- 
tions to virtual Z-boson production at the Tevatron, with Q = 910 GeV. The factorization scales 
indicated by the plateaux are 0.28Q for initial gluons and 0.53Q for real and virtual gluons. 
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FIG. 16: The factorization-scale dependence in the collinear scheme of the LO and NLO contri- 
butions to virtual Z-boson production at the LHC, with Q = 6.5 TeV. The factorization scales 
indicated by the plateaux are 0.28Q for initial gluons and 0A9Q for real and virtual gluons. 
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In Figs. n~3HT6l we show the factorization-scale dependence of the virtual Z-boson pro- 
duction cross section for higher values of Q, at both the Tevatron and the LHC. As Q/\fS 
increases, the factorization scales indicated by the plateaux decrease; this decrease is mod- 
erate for real and virtual gluons, and more pronounced for initial gluons. As the scale 
indicated by the initial-gluon plateau decreases, we observe a corresponding decrease in the 
scale at which the initial-gluon correction is minimized, which supports our argument for 
this scale choice. At large values of Q/y/S, the scales indicated by the initial-gluon plateau 
and the real- and virtual-gluon plateau differ by nearly a factor of 2. This is reflected in 
Figs. [15] and [HJ where the initial-gluon correction is small for /i pa Q/4, while the real- and 
virtual-gluon correction is small for p pa Q/2. When choosing a single factorization scale, it 
is necessary to compromise between these two values. Since the initial-gluon correction is 
dwarfed by the real- and virtual-gluon correction, for practical purposes one should choose 
the scale indicated by the real- and virtual-gluon plateau, /i ~ Q/2. 

In Figs. [17] and [18], we examine Drell-Yan production of a virtual photon in a fixed-target 
experiment with a beam energy of 800 GeV (yS = 38.8 GeV pp). Once again, the NLO 
corrections are small if the factorization scale is chosen according to the plateaux. Fig. [18] 
in particular illustrates the importance of this choice: if one instead chose \x = Q, the NLO 
correction due to real and virtual gluons would equal 35% of the LO cross section. 

IV. CONCLUSION 

In this paper we have introduced the collinear factorization scheme, and argued for a 
method to choose the factorization scale in that scheme. We applied this to the Drell- 
Yan process, and found that the explicit NLO corrections are small, less than 10%, for a 
wide variety of machine energies and values of the vector-boson invariant mass Q. The 
factorization scale indicated by our method is p pa Q/2 in the collinear scheme. This is not 
necessarily the correct scale in other factorization schemes, such as the popular MS scheme. 
We hope to adapt our method to choosing the factorization scale in the MS scheme. 

A complete calculation of a NLO cross section in the collinear scheme requires the LO 
partonic cross section to be convolved with NLO collinear-scheme parton distribution func- 
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FIG. 17: The factorization-scale dependence in the collinear scheme of the LO and NLO contri- 
butions to virtual photon production in pp collisions at y/S = 38.8 GeV, with Q = 5 GeV. The 
factorization scales indicated by the plateaux are 0.48Q for initial gluons and 0.64Q for real and 
virtual gluons. 
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FIG. 18: The factorization-scale dependence in the collinear scheme of the LO and NLO contri- 
butions to virtual photon production in pp collisions at y/S = 38.8 GeV, with Q = 20 GeV. The 
factorization scales indicated by the plateaux are 0.27Q for initial gluons and 0A8Q for real and 
virtual gluons. 
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tions. 4 Since these are not available at present, the numerical results presented in this paper 



have been computed using the LO CTEQ6L1 parton distribution functions [161 ] . The to- 
tal NLO correction then consists of the explicit radiative corrections plotted in Sections [UJ 
and IIH] together with the implicit correction resulting from the replacement of LO parton 
distribution functions with NLO parton distribution functions in the evaluation of the LO 
cross section. Without a NLO set of parton distribution functions in the collinear scheme, 
we cannot evaluate this implicit correction. However, if the factorization scale is chosen in 
such a way that the explicit correction is small, any large correction must be contained in 
the implicit correction. This suggests that a good approximation to the NLO cross section 
may be obtained by evaluating the LO cross section with NLO collinear-scheme parton dis- 
tribution functions, with the factorization scale determined as above. Such a set of NLO 
parton distribution functions could be generated once all the collinear-scheme counterterms 
are known. 
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